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Abstract 



One-loop QED corrections to the differential width of radiative muon decay are 
considered. Results can be used to analyze high statistics data of modern and future 
experiments. 
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1 Introduction 

Since the discovery of muon in 1937, the studies of its properties were always very 
important for the progress of the elementary particle physics. Nowadays, such precision 
observations like the muon life time and the muon anomalous magnetic moment are 
important for the checks of the Standard Model and searches for new physics. Besides 
many others, the process of radiative muon decay, 



is investigated in the modern experiments. In particular, the set of data from the 
PIBETA (tt/3) experiment [1] at the Paul-Scherrer Institute contains a considerable 
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amount of these decays. Accurate measurements of the process provide interesting in- 
formation about the structure of weak interactions. 

In this paper we construct an advanced theoretical prediction for the differential dis- 
tribution of process (1). Our calculations of radiative corrections (RC) allow to reduce 
the theoretical uncertainty. That makes it possible to perform precision comparisons 
with the experimental data and potentially look for new physics or rule out certain 
extensions of the Standard Model. 

In the limit of small energy loss (carried away by the neutrinos), radiative corrections 
to the process were considered in Ref. [2]. In this limit the standard decay produces a 
background to the searches for the neutrinoless decay fj, — > ej. 

In this paper we will consider the general kinematics assuming that the energies of the 
final state electron and photon are above of a certain threshold and the angle between 
their momenta is not small (see Sect. 3.3). The tree-level distribution and the notation 
are introduced in the next Section. Then we consider different RC contributions. In 
Conclusions we present some numerical results and estimate the theoretical uncertainty 
in description of the radiative muon decay. 



2 The tree-level distribution 



Within the Fermi model of four-fermion interaction, the differential width of radiative 
muon decay was first considered in Refs. [3,4]. Accurate formulae including the terms 
suppressed by the factor (m e /m M ) 2 were recently presented in Ref. [5]. We checked that 
their results coincide with the relevant contribution, which have appeared in calculations 
of exact one-loop radiative corrections to the muon decay spectrum [6]. At the Born 
level the differential distribution of the electrons and photons of the process (1) has the 
form 
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where Gp is the Fermi coupling constant; m e and m M are the electron and muon masses, 

— * 

respectively; P M is the muon polarization vector; x and y are the electron and photon 
energy fractions in the muon rest reference frame, x = 2E e /m^ and y = 2i? 7 /m /1 ; by 
p e and p 1 we denote the unit vectors in the directions of motion of the electron and 



photon, p e = p e /\p e \ and p 1 = p 7 /\p, 
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cos(p e p ). Functions F(x,y,d), G(x,y,d), 



and H(x, y, d) can be found in Appendix of Ref. [5]. 
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In what follows we will concentrate on the case of unpolarized muon decay, since it 
is the one measured in the PIBETA experiment. In the unpolarized case only three 
variables are relevant and the tree-level distribution can be represented as 



d 3 r 



Born 
unpol. 



Ot 



r Q —f3F(x,y,d). (3) 



da; dy dc 8ny 

Model independent parameterization of four-fermion interaction (see Particle Data 
Group [7]) leads to the appearance of two additional contributions. One of them is 
proportional to the difference (1 — 4p/3), which describes the deviation of the Michel 
parameter p from its value in the Standard Model. And the other one contains param- 
eter fj, which is a positive semi-definite quantity (see Ref. [8]) 

V= (\9 V rl\ 2 + \9lnl 2 ) + l(\9 S LR + 2g T LR f + \g s RL + 2g^) + 2(|^| 2 + |^ L | 2 ), (4) 

where g^LR are the right-left (RL) and left-right (LR) coupling constants, which pa- 
rameterize non-standard scalar (S), vector (V) and tensor (T) four-fermion interactions. 
In principle, one can look also for other exotic interactions, e.g., for the ones mediated by 
antisymmetric tensor fields [9]. Extraction of fj from the experimental data potentially 
can put strict limits on physics beyond the Standard Model. 



3 Radiative corrections 



New precision experiments call for an adequate level of accuracy in theoretical predic- 
tions within the Standard Model. Effects of higher orders of the perturbation theory 
become important. Here we will consider the first order QED radiative corrections. As 
usually, we separate them into three parts: i) emission of an additional soft photon; ii) 
effect due to one-loop virtual photonic correction; Hi) emission of an additional hard 
photon. Note that all the relevant pure week corrections (like loop insertions into the 
W^-propagator) are included into the Gp coupling constant [10,11], which is measured 
directly from the muon lifetime. Effects of strong interactions in the process under con- 
sideration are negligible for the moment. They start to appear only at the order O (a 2 ) 
through hadronic vacuum polarization. 



3. 1 Soft Photon Contribution 



We assume, that emission of an additional soft photon of energy below certain threshold 
is not distinguished by the experiment from the tree-level process (1). The energy of 
the soft photon, uj 2 , is limited by the parameter A: 
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The corresponding correction can be factorized out in front of the tree-level differential 
distribution: 
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where A is a fictitious photon mass; L is the so-called large logarithm, L = ln(m 2 /m 2 ) 
10.66 ; the dilogarithm and the Riemann zeta-function are defined as usual: 
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Quantity 5 Soft coincides with the corresponding factor, arising in the correction to the 
non-radiative muon decay (see e.g. Ref. [6]). Expression (6) takes into account the de- 
pendence on the electron mass exactly. Omitting small terms proportional to (m e /m^) 2 , 
we get 
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3.2 One-loop virtual correction 



Here we will consider the effect of one-loop photonic corrections. Some representatives 
of the relevant Feynman diagrams are given in Fig. 1. There are two diagrams of class 
(a) with photon emission from an external leg (electron or muon line). In the same 
way the two box-type diagrams of class (6) describe real photon emission from virtual 
electron and muon propagators. Diagrams of classes (c) and (d) give corrections to 
photon radiation from a single leg. To get the corresponding correction to the muon 
decay spectrum we have to multiply the complete set of amplitudes of classes (a — d) 
by two tree-level amplitudes, describing single photon emission. In our calculations we 
followed the procedure which has been applied in Ref. [12]. 
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(d) 




Fig. 1. Types of Feynman diagrams for radiative muon decay with one-loop RC. 

The standard technique for one-loop integration was used. The list of relevant integrals 
is given in Appendix A. To eliminate the ultraviolet divergences we applied renormal- 
ization of the masses and wave functions of the electron and muon. Note that this is 
enough in the case of muon decay (see Ref. [13,14]), contrary to the general case of 
the Fermi four-fermion interaction. An analytical result for the virtual correction was 
obtained. We do not give the full formula here, since it is rather long. 



3. 3 Emission of an additional collinear hard photon 



Events with registration of two hard photons are supposed to be rejected by the ex- 
perimental event selection. But if the additional photon is emitted at a small angle 
with respect to the momentum of the outgoing electron (positron), the former is not 
recognized by a calorimetric detector as an independent particle (this can happen if 
there is no any considerable magnetic field in the detector volume). So, for the so- 
called collinear photon emission, one observes an effective electron with the energy 
and momentum composed by the sum of the corresponding quantities of the photon 
and the bare electron. Let us assume that this kind of calorimetric registration hap- 
pens in the experiment, if the angle between the electron and photon momenta does 
not exceed a certain value 9 , which plays the role of a small parameter. We demand 
m e /m^ < < 1- Typical experimental values for this parameter, a few degrees, sat- 
isfy our conditions. On the other hand, the angle between the observed photon and the 
electron should satisfy the condition 6 = p e p^ > 9 . 

According to the general factorization procedure, we can represent the result for the 
contribution of collinear photon radiation as the product of two factors: 



j3-pH— coll j3-pBorn 
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Rcoll = ^ J y{[l + (l-^) 2 ](^ + 21nx-l + ln| + 21n(l-^))+^}. 

A/x 

The tree-level radiative muon decay (with photon emission at large angles with respect 
to the electron momentum) serves as a short-wave sub-process. Emission of a collinear 
photon by the outgoing electron serves as a long-wave sub-processes. The formula for 
the collinear radiation factor agrees with the one in Ref. [15]. 

Integration over the energy fraction of the collinear photon, z, gives 



R. 



coll 



a 
2^ 



o2 

L + 2\nx- 1 + ln^ 
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Note that the lower limit of the collinear hard photon energy fraction is adjusted to the 
upper limit of soft photon emission. 



4 Results and conclusions 



Summing up the contributions of soft, virtual, and hard collinear photonic corrections 
we receive the final answer for the first order radiative correction to the process (1). 
Here is our result for the corrected distribution, which substitutes function F(x, y, d) 
from Eq.(3): 



F Corr -(x, y, d) = F{x, y, d) (l + ^A(x, y, d)) + ^B F (x, y, d), 

A{x,y,d) =2 ln^ln a: -In A) - 2 In A - ^ In ^ + ^ln ^ - 2 hix^ . (11) 

We presented explicitly only the factorized part of the correction. The remaining non- 
factorizable part, Bp(x,y,d), is rather long. We use it in a FORTRAN code for numer- 
ical estimates. Expressions for the radiatively corrected functions G Corr '(x,y,d) and 
H Corr '(x,y,d) have exactly the same form as Eq. (11) with the trivial substitutions: 
F — > G(H) and B F (x,y,d) — > B G ( H )(x,y,d). The most important factorized part of 
the correction, A(x,y,d), is universal for all the three functions. 

It is worth to note that all the leading logarithm terms were factorized in each of 
the contributions, but they cancel out in the sum in accord with the Kinoshita-Lee- 
Nauenberg theorem [16,17]. Moreover, all the dependence on the parameters A and 6q 
is contained in A(x,y,d). 

In Fig. 2 we plotted the Born-level differential branching ratio of the radiative muon 
decay for a fixed value of c, 
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Fig. 2. Differential branching ratio versus electron energy fraction y for three different x- values 
with fixed c = 0.5, A = 0.01, = 3°. 

The relative contribution of radiative corrections is illustrated by Fig. 3, 



F c ™ix,y,d)-F(x,y,d) 
F(x,y,d) 



100%. 



(13) 



The dependence on c value of <5 RC is rather weak. For given values of c and x the 



x-0.2 
x-0.4 
x=0.6 



O 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

y 

Fig. 3. Relative contribution of radiative corrections versus the electron energy fraction; pa- 
rameters are the same as in Fig. 2. 

maximal value of y is defined by the kinematics: 



|/max 



1 — x + ml/m 2 ^ 
l-x(l-Pc)/2 ' 



(14) 
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For the given set of parameters, the factorized part of the correction dominates and 
gives about 4/5 of the total effect. 

To illustrate also the case of 100% polarized muon decay we present in Fig. 4 a plot 
for the relative contribution of radiative corrections for a set of fixed variables. Namely, 

P/j,P e = 30°, P (i V 1 — 60°; c, 9 and A are the same as in Fig. 2. Quantity 5^. is defined in 
analogy to Eq. (13) by adding the relevant contributions of G and H functions according 
to Eq. (2). 



So -4.8 



x o 2 

x-0.4 
x=0.6 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



Fig. 4. Relative contribution of radiative corrections for the case of polarized muon decay 
versus the electron energy fraction. 



Thus we presented the calculation of one-loop QED corrections to the differential dis- 
tribution of unpolarized muon decay. Our FORTRAN code is available upon request from 
the authors. The results can be applied also for the decays r — > /iz/ M z/ T 7 and r — > ev e v T ^. 
The theoretical uncertainty of the spectrum description is defined by higher order QED 
radiative corrections (EW and QCD effects are negligible compared to the QED ones). 
As a rough upper estimate we can consider the relative contribution of the omitted 
higher order terms to be about (5 RC ) 2 < 3 • 10 -3 , which is small compared to the present 
experimental precision. In principle, one can easily get the most important higher order 
terms with logarithms of A and 6 by means of soft and collinear approximations. 
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Appendix A. 
One— loop integrals 

Here we give the list of integrals over the loop momentum k%, which are used for 
calculations of the virtual loop contribution. The notation for the loop integrals is as 
follows: 




(A.l) 



where (i — I) are the denominators of the relevant propagators: 
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(0) = k\ - A 2 , (1) = k\- 2k igi , (2) =k\ — 2k lP , 

(3) = k\ - 2h( qi + k 2 ) + X e, (4) =k\ — 2h(p - k 2 ) - Xl » 

Xe = 2k 2 qi, X» = 2k 2 p- (A.2) 

Tensor and vector integrals are decomposed as 

i» v — n ^j9 i^^tpp i//r« j^h^k 1 "!^ +f// + nVU M 

1 ijkl — 9 1 ijkl + P P I ijkl + Ql ( ll I ijkl + K 2 K 2 2 ijkl + \P Hl+P Ql) I ijkl 

+ (p"fc* + + (« + t&Wlgu, 

Iijki = P^Iijki + Qi lijki + k% lijki- (A-3) 

We need different integrals with 4, 3, and 2 propagators. In the integrals below we 
dropped the dimension and put m M = 1 (the dimension can be restored by multiplying 
by the required power of m M ). In some cases we express tensor and vector integrals 
through a combination of more simple ones. These relations can be obtained by mul- 
tiplying the tensor and vector integrals by certain particle momenta and subsequent 
cancellation of denominators, where possible. 

Ultraviolet and infrared divergences are regularized by introduction of a cut-off, A, and 
an auxiliary photon mass, A. They appear in the integrals in two logarithms: 



La — In —r, A>m^ (A.4) 
ml 

L\ = In —^i A<m e . (A.5) 
m\ 



The short notation z — \xyd will be used below. The relevant tensor integrals are 
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y — z\ 1 — x 



The following vector integrals were used in our calculations: 



-^0124 — ~ (~ 2-^0124 ~~ X ^0124 + ^124 ~ ^014^ , ^0124 ~ ~ 2/^0124 + Uhl2A + I{ 

- In 2 x + ln(x — z + y) ^ln(x — z) — ln(l — x + z — y)j 



1/1 



l 0124 — o 2 



+ lny(hix + ln(l - y) - ln(x - z) - ln(x - z + y)^ + ^ln 2 y~C(2) 
/y + x — 



+Li 2 



— Li 9 



z - y 



x 



- Li 2 (x - z + y) + Li 2 (y) - Li 2 (1 - x) , 



•^0123 ~~ ~ (^0123 ~~ ~ \Jo23 ~ h:2:>, \ ~ - h)V2:', + h)V2 ~ h)23 J - 
^0123 — ~ ( — 2-^0123 ~~ -^023 + h23 ) 



13 



-'0123 



' z ( z - X y) i^2 X y Zl ^ ~ \ X Vih\2 - J 023) ~ 7^(4l3 ~ ^123 ) 



+z(Iqi2 - -^023) - ^0123 + 2^ 2 (^023 - ^123)^, 

to 1 ( T n In x \ riD 1 /. lux \ 

J oi2 = [L + ]nx + - , ij 12 = lnx-- , 

XV 1 — xj XV 1— X/ 



'013 



2-^013 — zlmz + Im — I 



y 



-^014 — 

(x — z) 2 



013 ~r J01 — J 13 

x — z 



I? 



013 



"Z/oi3 + -^01 



'03 



/ X — Z \ 

L ^ ln(x + y - z) + In yj 



+ ln(x + y-z)(2- 



y 

z — x + 1 



+ 



1 



y y(l-x-y + z) 1-x-y + z, 

— 2 lny + 2 my ln(x — z) — \i\y ln(l — y) + 2 lny ln(x + y — z) — ln 2 (x + y 
—2 ln(x + y — z) ln(x — z) — ln 2 (y) + ln(x + y — z) ln(l — x — y + z) 

' x + y - z" 



-2Li 9 



y 



-^014 



x — z 



f + Li 2 (x + y-z)-Li 2 (y) + 2C(2) 

y , x+y-z 

■my — ; — m(x + y — z) 



1-y 1—x—y+z 
J 023 = ( J 023(^ + 2/) + ^os - ^23) , 

^03 + y) - / 23 (x + y - 2) + 2z/ 023 - 2/ 02 ) , 



-'023 



(x + y) 2 - 4z 



/ P 24- ~ 

j k __1 

024 — 



-^123 — 



lny-Li 2 (l-y) + C(2) 



2 - lny 1 + 



1 



i-y 

x + y-z 



+ - Li 2 (l-y)-C(2) 



r 1 ? 

J 123 



y- z 
1 

y - 2; 
x + y — z 



ln(x + y — 2;) + 



x 



1 -x 



lnx 



1 — x — y + z 

L^— lnx + ln(x + y — z)^j + ln 2 (x + y — z) — lnx — 

ln(x + y — z) — ln(x + y — z) ln(l — x — y + z) 



In 



1 — x — y + z 
-Li 2 (1 - x) - Li 2 (x + y - z) + C(2) 
1 



r - 

J 123 ~~ 



+ 



y-z 
x 



y-z 



2 + L[ -1 — lnx H — ln(x + y - z) ) 

V y-z y-z / 

^C(2) - ln(x + y - z) ln(l - x - y + z) + ln 2 (x - z + y) 



-Li 2 (x + y — z) — Li 2 (1 — x) + 2 lnx — ln 2 (x) 
2x 1 



- ln(x - z + y) 1 + 



y — z 1 — x — y + z 
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J 124 



y-z 



Inx + : X ■ — ln(x + y — z) + Li 2 (1 — x) 



l-x 



1 — x — y + z 



-Li 2 (1 - x - y + z) 



T q - 

J 124 — 



1 



'124 



y-z 
1 



x 



l-x 



Inx 



y-z 



2 + 



x 



Inx 



x + y - z 
1 — x — y + z 
x + y - z 



1 — x 



X 



y-z y-z 
x + y — z 



ln(x + y — z) 
ln(x + y — z) 



Inx — 



1 — x — y + z 



y — z VI — x 
2 — x i 

' Li 2 (1 - x - y + z) - Li 2 (1 - x) 



ln(x + y — z) 



J oi 



-^04 



y-z 
1 1 



1 



-- + - [L A + L , J£ 2 = -- + -L A , J 9 3 = - + - L A + L-ln^ , 



1 1 / r r 

1 1 

- + - L A - Inx + 



lny 



4 ' 2 
1^)' 



(l-x 



1 . 1 

mx + 



1-xJ' 



T p - 



r q 



X 



1 — x (1 — x)'- 



lnx^j , 



L A + L-^-, I k 13 = -(L k + L-°- , II 



'-14 — J 14 



J 14 r>> 



-^14 — -^235 -^23 _ -^23 _ -^23 ~~ 7T' -^24 -^A ~ ; -f 2 4 — ~ ~ ( L 



J A — - 



27' 



-^23 



1 1 

4 + 2 



1 — x — y + z 



+ ln(x — z + y) ( — 1 + 



(1 - x - y + zf 



The scalar integrals read 



-^0123 — 



xz 



1 1 

-L\L — L\ Inx + -L 2 — L ln(x + y — z) 



+Z/lnx + L\wz + 2 Inx In z — ln 2 (x — z + y) — 2Li 2 



-^0124 — — 
xy 



-L\L — L\ Inx L 2 — L ln(x — z + y) 



+L In y + 2 In a; In y — In 2 (x — z + y) — 2Li 2 



z - y 



a; 



g - y 

X 



C(2) 



-C(2) 





1 


1 


-^012 — 








X 


.2 




1 


1 


^013 — 










.2 


-^014 = 







; LL A + lnxL A + -L 2 - In 2 x - Li 2 (1 - x) 



; (L + lnz) 2 -C(2) 



x — z 



L(ln(x + y — z) — lny) + ln 2 (x + y — z) + lnyln(l — y) 



+ ln(x — z) ln(x + y — z) — In y ln(x + y — z) — \ny ln(x — z) 
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+Li 2 

-^023 



x + y - z 



2Li 2 



+ Li 2 (1 - x — y + z) — Li 2 



z — x 



+ Li 2 (y)-2C(2) 



r(x + y) 



1 + r , 1 + r 

x + y) 



2 v " ' " y l-r\ 
1 . o 1 + r 1 , o 1 - / 
2"' 2 2 



2Li 2 1 



-Li 2 (1 - x - j/ + z) + ^ In 2 — ; In 2 - In ^-r- In 



(1 - r)(x + y) t 
1+r, 1 — r 



-2\n(x + y) In ln 2 (x + y) 



r = Jl 



Az 



(x + yf 



1 



'024 



'123 



Li 2 (l-y)-C(2) 
1 



y-z 



L(\n(x + y — z) — In x) + ((2) — Li 2 (x + y — z) — Li 2 (1 — x) 



In 2 x — \n(x + y — z) ln(l — x — y + z) + ln 2 (x + y — z) 
1 



'124 



y - z 



Li 2 (1 - x) - Li 2 (1 - x - y + z) 



ioi = L A + L + 1, J 02 = L A + 1, J 3 = + ^ + 1 - In x, 

Io4 = L A + 1 + In y, I12 = L A + 1 + In x, I 13 = L A + L - 1, 

1 — y 1 — x 

x ~\~ y — z 

hi = I23 = L A + 1 + ■ — ln(x + y - z), J 24 = L A -1. 

1 — x — y + z 



16 



